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In a previous paper [J. Phys. Chem. Solids 27, 975 (1966)] the authors have discussed symmetry 
conditions for equivalence of jump rates operative in point-defect motion between crystallographically 
equivalent general sites neighboring a trapping center. The treatment is here extended to equivalent 
special sites to show that the total number of inequivalent jump rates can be expressed as (Nt/Nd)~~Qe~ 1 
where N t and N ( t are the orders of the symmetry groups of the trap and defect respectively and where 
q,. is the number of independently effective pairs of nonself-inverse symmetry operators (paired with 
their inverses) all in the symmetry group of the trap. The number q e can be obtained by counting any 
nonself-inverse operator pair, an element of which, when multiplied by any element of the defect 
symmetry, equals neither its own inverse nor an element of a pair previously counted. 

For thermally activated relaxation processes expressions apply for maximum and minimum 
numbers of jump frequencies involved in complete and partial relaxation processes. For complete 
relaxation the maximum number is the number of prime factors in the ratio of the order of the symmetry 
group of the trap to that of the defect group. The minimum number is the minimum number of gen- 
erators that will raise the defect position symmetry to that of the trap. 
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1 . Introduction 

In a recent publication [1] l we have discussed sym- 
metry conditions for equivalence of jump rates opera- 
tive in point-defect motion between crystallographically 
equivalent sites with a view toward an assessment of 
the extent to which such motion can account for relax- 
ation and diffusion phenomena observable in crystals. 
Closely related problems have been treated by A. S. 
Nowick [2] and in papers referred to by him. Our treat- 
ment does not depend on linearity of response to an 
outside influence and is believed to have rather general 
applicability. In the quoted paper [1] it was applied to 
a highly specialized situation governed by four restric- 
tions; the partial removal of the second restriction 
being the contribution in this paper. 

Restriction 1: The point defect, though itself rela- 
tively mobile, was held to an immobile trapping center, 
the defect remaining free only to jump between sym- 
metry-equivalent sites on a shell neighboring the 
trapping center. The effect of this restriction was to 
reduce a problem concerned with an infinite number of 
possible defect sites to one concerned with the finite 
point groups of the trap, G t and G\ in the initial and 



1 Figures in brackets indicate literature references at the end of this paper. 



strained conditions, respectively (all basic notation is 
summarized in table 1). The corresponding point-group 
orders are designated N ( and N' t . Both groups are sub- 
groups of the corresponding crystal space groups, and 
of the crystal point groups G and G' of order N and N' , 
respectively. 

Restriction 2: The point group of the point-defect 
site in the unstrained crystal, Gd, was restricted to the 
identity 1. In the strained crystal the point group of 
the defect, G^, was necessarily 1, too, because by 
simple symmetry postulate [1] G' d QGd. We here relax 
this restriction to allow the site symmetry of the defect 
to differ from the identity and we refer to its symmetry 
as symmetry of specialization. Gd, of order 7V r /, must 
be a noncentrosymmetric subgroup of G t because the 
trapping center neighbors the defect on only one side. 
By the correspondence of the shell of sites around the 
trap with the construction sphere for a stereographic 
projection, we deduce that the only permitted Ga 
are the ten symmetries achievable on the surface of 
a sphere (see tables 3 and 4 in reference [3]), which 
in turn are isomorphous with the ten two-dimensional 
diperiodic point groups [4]. Although in three dimen- 
sions the presence of the trap imposed the limitation 
of lack of a symmetry center, this limitation does not 
apply to the isomorphous two-dimensional groups. 
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Table 1. Symbols used for point groups, their orders as well as for 
individual operators, defect sites and paths 



Crystal 


Point 
group 


Order 


Symmetry 
operator 
in group 


Crystal 

Trap 

Defect 


Unstrained 
Strained 

Unstrained 
Strained 

Unstrained 
Strained 


G 

C 

G t 

G d 
G' d 


N 

N' 

N t 

N' t 

N d 


S 
S' 

s t 

s d 
s d 



G d CG d GG t QGDG' D G[ D G d and G t D G\ 

where D indicates a subgroup relationship of the right group to the 
one written to the left of that_symbol._ 

An S could be a 1, 2, 3, 4, 6, 1, 2, 3, 4 or 6 axial operator to the first 
positive or negative power. 

S~ x is specifically an inverse of an S. 

S is specifically an inversion axis operator. 

S is specifically a rotation axis operator. 

\S\ is the order of the axis. 

Additional suffixes 1, 2, 3 . . .are used to distinguish different 

operators. 
[uvw] is a direction e(uvw) where e means "an element of." 
S[uvw] is the direction e{uvw) to which [uvw] is transformed by 

operator S. 
S2Si[uvw] =S 3 [ai;M;]when the successive operations Si followed by 

S 2 is equivalent to the single operation of S ; $. 
[uvw] —> S[uvw] is a denned path from a site defined by [uvw] to one 

defined by S[uvw]. 
[uvw] A S is the angle between [uvw] and the direction of S. 



Indeed the two-dimensional environment of a defect 
on the surface of the shell and on its tangent plane may 
be two-dimensionally centrosymmetric. 

There exist alternative ways of expressing all pos- 
sible symmetry operators S. Most systems commonly 
used allow, in addition to 1, any of the crystallo- 
graphically permitted rotation axes 2, 3, 4, or 6. To 
complete a self-consistent set of distinct possible 
operators S in G we will use inversion axes 1, 2, 3, 4, 
and 6. Where the order of the operator, \S\, defined by 
the smallest positive integer for which S' 5 I = 1, is 
greater than 2 (that is for 3, 4, 6, 3, 4, and 6) the operator 
is nonself-inverse and its inverse is distinct. Never- 
theless, if S is an element of G, S~ l is, too, i.e., 
SeG 4 ^* S^eG. Of the inversion axes only 2, equivalent 
to a mirror plane, m, provides a possible SdtGd. 2 differs 

from all rotation-axial Sd, denoted S<j, in that 2d pro- 
vides a locus of an infinite number of coplanar direc- 
tions corresponding to defect sites all of which have 
Gd 3 2, whereas a crystal can only possess a small 
number of sites with rotational symmetry of specializa- 
tion. (It is not implied that a defect with 2d has a degree 
of freedomof motion alongthe mirror plane perpen- 
dicular to 2d. The sites on 2d are not in general sym- 
metry related and so can have differing free-energy 
characteristics.) The possible G& are then the five 
Sd alone plus the same five with a perpendicular 



2d axis (Sd/2d). Thus only ten defect groups can occur: 
1, 2, 3, 4, 6, m, 2m, 3m, 4m, and 6m. The effect of the 
second restriction imposed in our previous paper [1] 
was further to reduce consideration to equivalent 
general sites of the point defect. With G d =l, neces- 
sarily, the number of equivalent defect sites was always 
N t . 

Each site can be defined by a [uvw]e(uvw) where the 
u, v, w numbers need not be simple integers. Equality 
of two directions represented by [u\, V\, W\\ and 
[u 2 , v 2 , w 2 ] is assured by the conditions that the magni- 
tude and sign of any two of the three ratios -, — , — be 

v w u 

the same in the two cases and that the sign of one of 
the three symbols u, v, w be the same. Every element of 
(uvw) can be described as St[uvw], the direction gen- 
erated by S t eGt acting on [uvw]. A path indicated by 
[uvw] — > St[uvw\ will refer to a path from a defect site 
[uvw] to another site St[uvw]. Of the infinite number of 
such geometric paths a choice can be made to fulfill 
some specific conditions. If the notation were applied 
to crystal-chemical networks, it might be adequate to 
designate the straight-line path. In other applications 
a great-circle path is satisfactory. In our work here, we 
mean the path which has the largest jump rate, which, 
in turn, approximates that of minimum free-energy 
barrier. The path for largest jump rate need not be 
straight, nor lie on the surface of the shell containing 
the sites and has no single geometric representation 
to justify fully the use of the singular word "path" 
which, nevertheless, we use. One may think of a path 
with a finite not-sharply delineated width. The removal 
of restriction 2 will necessitate a further reexamination 
of the concept of the path (see sec. 2). 

In our previous work [1], for which always Gd=l 
and for which the number of sites, consequently, 
equaled the number of elements in (uvw), i.e., TV*, 
the number of nonzero paths emanating from [uvw] 
was always N t — 1. Symmetry requires that these 
paths from [uvw] are fully representative of all paths 
relevant to (uvw). This rule will apply equally for 
N d >l. 

The aim of the present paper is then to remove 
this second restriction, but we simplify our task by 
restricting discussion to noncubic trap sites G t . 

Restriction 3: We have assumed that (in our new 
nomenclature) 

[uvw] — » St[uvw] = [uvw] —> S y l [uvw] (1) 

even when \St\ > 2 that is when 

[uvw] —> St[uvw] ^ [uvw] —» Sj l [uvw]. 

This eq (1) will hold for most physical models one can 
invent. In our previous paper [1] we have in fact shown 
that the equation must hold for all steady-state con- 
ditions under which the principle of detailed popula- 
tion balance is applicable. Our treatment, though in 
some ways complicated by this restriction, would 
lose likely contact with possible experiments if we 
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removed it. We, therefore, retain the condition im- 
posed by eq (1) throughout this paper. 

In our previous paper fl] we found that for each 
pair of nonself-inverse operators S and S _1 , the num- 
ber J of nonequivalent paths from [uvw] was reduced 
by one. If q is the number of such S pairs (\S\ > 2) 



J = N t -q-\ 



(2) 



Even after completely upsetting the population 
in equivalent sites by an externally imposed influence, 
relaxation by a thermally activated process does not, 
in general, require all the nonequivalent jump rates 
to be activated. It is in general impossible to state a 
number of needed jump rates even for a given G t . It 
was, however, possible [1] to derive simple expres- 
sions for the minimum number J min , and the maximum 
number ./ max needed completely to relax a system. 
The minimum number, J mxn , of nonequivalent jump 
rates needed completely to relax a system in which 
the population of defects over equivalent sites had 
been completely upset was the rank of G t (the mini- 
mum number of operators whose products and powers 
constitute that group). In most point groups the ex- 
perimentally observable number of jump rates could 
be larger than y min depending on the relative magni- 
tude of the jump rates. There was, however, a usefully 
restrictive maximum number / max of nonequivalent 
jump rates equal to the number of prime factors in 
N(. Similarly, simple rules could be formulated for 
partially upset populations. If N' t of the N t initially 
equivalent sites remained equivalent through a postu- 
lated sequence of influences, the minimum number of 
jump rates, J min , needed for relaxation, was the mini- 
mum number of nonequivalent symmetry operators 
that would raise G\ of order N\ to the initial point group 
G t . 

The maximum number was the number of prime 
factors in (Nt/N f t ). We presented [1] tabulations of 
maximum and minimum numbers of jump rates. 
Experiments designed to observe these points must 
be carried out at temperatures high enough so that an 
adequate number of jump rates are above the level 
of observation, so that point defects can attain equilib- 
rium distributions in finite time. 

Restriction 4: Along paths emanating from a given 
[uvw] predictable cross-over positions may occur with 
equivalent paths emanating from other elements of 
(uvw). Some crossings must be exact, some may only 
be approximate, because of limited variability of the 
radial distance from the trap. As a result some paths 
must have equal, approximately equal, or lower free- 
energy barriers than others also emanating from [wi;m;]. 
These inequalities have not been discussed in detail, 
nor are they further discussed in this paper. It should 
be realized, however, that the effect of these inequali- 
ties will be to lower the number of experimentally 
observable distinct jump rates. Accidental equalities 
of jump rates under specific conditions and within any 
given experimental precision need not be discussed 
in the present context. 



2. Consequences of Symmetry of Specializa- 
tion for the Defect Sites 

By analogy with the definition of general posi- 
tion in crystallography as the set of symmetry-related 
positions each of which is left invariant by no other 
symmetry operator S than the identity, 1, we will use 
"general direction" to define the (uvw) set if a [uvw] 
e (uvw) is not invariant under any operator St other 
than the identity. Similarly "special direction" will 
signify the (uvw) set if at least one St (other than the 
identity, 1) will leave a given [uvw] invariant. The 
special direction characterizes the symmetry prop- 
erties of a special position set of a point defect con- 
fined to a trap. The stereographic representation of 
a special direction set (uvw) can easily be visual- 
ized with the help of the steoregraphic representation 
of the general direction set for the given G ( [3]. The 
positions of the general set (uvw) are moved, in con- 
formity with the symmetry requirements, on the 
stereogram into coalescence on the stereographic 
representation of the defect symmetry group G<t. 

(tt — 1) paths will then emanate from any [uttt#]. 

The direction S d [uvw] = [uvw] so that generally 
for SttSd'. 



S t [uvw] = S t S ( i\uvw\ 



(3) 



where the product S t S ( i signifies the operation re- 
quired by Sa followed by that of St, a product which 
must equal an operator SficG*; where S n ^ St unless 
Sd=l. We, therefore, note that for special positions 
any one direction of the set (uvw) can be written in 
terms of a given [uvw] in Nd ways. All will be of the 
type St[uvw]. A given path [uvw] — > S ( [uvw] can also 
be written in Nd ways referring to identically the 
same path by changing only the operator S,._This is 
illustrated for the case G, = 6w2, G,/ = 2 1 (2 ± signi- 
fies 2 ±6) in figure 1 with its accompanying table. 
The concept of the path itself needs to be examined 
once more also by reference to figure 1. Reference js 
made to the path [uvw]^_2\\[uvw] for which S d = 2 ± 
leaves both [uvw] and 2||[z/i;w] invariant. For all 
models for which the geometric path is permitted to 
depart from the mirror plane of specialization, there 
will be symmetrically equal geometric realisations 
of that path on the two sides of the symmetry plane, 
2d. The 2[io!oi axis relating [aittc;] to 2\\ [uvw\ has a sim- 
ilar effect on any path not going through the angle 
bisector. Although we are aware of this plurality 
of some geometric paths [uvw]—*St[uvw] we feel 
justified in continuing to use the singular word path 
giving rise to one jump rate. For a given free-energy 
barrier a geometric path plurality will slightly increase 
the jump rate, but that enhanced jump rate will 
be a characteristic of the system, namely, that ap- 
propriate for the path [uvw\ — > St [uvw\. The geometric 
plurality of paths does not introduce complexities 
any more than did the path width discussed in the 
introduction. In fact a wider path with energetically 
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FIGURE 1. Equivalence of paths from a special position on a mirror plane in point group 6m2, 
for a defect on a mirror plane, m = 2. 



234 



shallower sides will in the same way lead to lower 
jump rates. 

In the next section we shall obtain the number 
of independent paths by a procedure involving the 
counting of the effective number of nonself-inverse 
operator pairs. We conclude the present section by de- 
scribing an alternate procedure involving the counting 
of certain distinct operator product sets; this alternate 
procedure is perhaps not as easy to apply to specific 
cases but provides a different point of view. We begin 
by noting that we need only consider the paths emanat- 
ing from any one [uvw] and that all paths emanating 
from a given [uvw] are of the form [uvw] — * S t[uvw] 
where S t # Gd because [uvw] is invariant under all 
operators in Ga> For two paths to be equal as ex- 
pressed by [uvw] -» S n [uvw]= [uvw] — » S t2 [uvw] the 
operators must satisfy either of two conditions. First, 
equality of paths is assured if S^ = S^ 1 and Sti is non- 
self-inverse; this is the only case which can occur for 
general [wt>w;] and was previously treated for general 
sites. Second, if SdtGd and 



S t 2 = SdSn or S/2 — SdS n l 



(4) 



equality of paths is assured (note that the first con- 
dition is contained in the second by taking Sd to be 
the identity). This second condition holds because 
an operator Sd which satisfies this condition will 
leave [uvw] unmoved and will take S n [uvw] into 
Sn[uvw]. By considering each of the possible opera- 
tions Sd in turn one can see that each point of the 
path [uvw] —> Sn[uvw] is carried into the corresponding 
path [uvw] — * StJ[uvw] so that the paths are equivalent 
and must have equal jump rates. This argument 
justifies the following procedure for determining 
the number of independent paths, J: list all of the 
StiGd. For each St form the set of operator products 
SdSt and SdS^ 1 using all SdtGd- J is equal to the number 
of distinct sets. 



3. The Number of Distinct Jump Rates in 
Terms of the Effective Number of Non- 
self-inverse Operators 

In this section we shall first show that only non- 
self-inverse operators need be considered in deter- 
mining the number of equivalent paths because 
even though self-inverse operators do in certain in- 
stances relate pairs of paths emanating from special 
positions their effect is always duplicated by the ac- 
tion of nonself-inverse operators. We then consider 
the fact that distinct pairs of nonself-inverse opera- 
tors sometimes coincide in their effect on paths 
emanating from special positions or fail to relate 
distinct paths so that the effective number of nonself- 
inverse operator pairs is sometimes less than the total 
number of distinct pairs. 

Two physically distinct paths emanating from [uvw] 
may be equalized by self-inverse elements of Gd 
acting directly on the paths or by the demands of 



nonself-inverse symmetry operators in G t (compare 
restriction 3 in the introduction). We can easily 
cite examples of the second kind of cause without 
the first being applicable (e.g., G t = 4/m; Gd=2; 
[uvO] -» 4>- l [uv0] = [uvO] -» 4[ut;0]). 

The two paths related by nonself-inverse operators 
(eq. (1)) need not be related by any element of Gd. 
Conversely, however, we can prove that two paths 
related by a self-inverse element of Gd are necessarily 
also related by a nonself-inverse operator as in eq. 
(1). Let Sd\ be a self-inverse operator of Gd which 
leaves [ai;w;] invariant; the path [am>] — ■> St[uvw] 
must be equivalent to the path [uvw] — » Sd\St[uvw]. 
If S ( [uvw] is also on Sdi then two paths connect 
the same directions and so are considered together 
as a single path as previously discussed and no 
equivalence of different paths results. If S t [uvw] 
is not on Sd then we can show that either S t is non- 
self-inverse or Sd\St is nonself-inverse. To prove this 
statement the Jour possibilities of self -inverse oper- 
ators 1, 2, 1, 2 for Sdi must be considered in turn. 
The possibilities 1 and 1 can be dismissed at once 
because the former has no effect and the latter cannot 
be an element of any Gd. For Sd\ = 2 any mirror 
plane whose normal is 1 1 to Sdi will leave [uvw] in- 
variant; a center, mirror plane, or other 2 axis _L Sdi 
will generate [Uvw\ still on Sdi. Thus the only self- 
inverse operators St which will move S t [uvw\ off 
Sdi are 2 or 2 axes inclined to Sdi whose product 
with Sdi is a nonself-inverse operator._ An exactly 
analogous argument applies for Sd2 = 2. Note also 
that if any \Sa\ > 2 for noncubic G t all elements 
of {uvw) are coincident with Sd. No SdtGd can then 
repeat any path [uvw] — * S f [uvw] other than into its 
inverse (i.e., the same path with the direction of 
motion reversed). We are, however, only concerned 
with the number of paths emanating from [uvw\. 
We conclude that any pair of paths related by a self- 
inverse operator is also related by a pair of nonself- 
inverse operators as in equation 1 and so turn our 
attention to examining the latter. 
If G t contains no \St\ > 2, i.e., when 4 = 0, no dis- 
tinct paths can be equal, so 7 = ^ — 1. The total 

d 
number of nonself-inverse operator pairs may not 

be effective in relating distinct paths when the (uvw) 
is a special direction set. For a general direction set 
(uvw): S t \[uvw] ^ S&luvw], However, for a special di- 
rection (uvw) it could be that S n [uvw] = S tl [uvw]. 
Two cases arise: (1) S n =S^ and S 2 tl [uvw]=\. This 
will apply if and only if [uvw] \\ S t (\S t \ >2). For 
noncubic G t all S, (|S,|>2) must be parallel to each 
other, and for (uvw) \\S t (\St\ >2) the number q e of 
independently effective pairs of nonself-inverse opera- 
tors is always zero. (2) The other special case is that 
in which S n and S t 2 are not a mutually inverse pair, 
yet generate the same path emanating from [uvw]. 
Equation 3, however, expresses all Nd possible ways 
of expressing a given path. It follows that the con- 
dition for this special case is: 

Sti — StzSd- 
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However, for Sn||S*2, all 2 and 2 J_ S, relate all S t 
to Sf l only. For 2||S/, i.e., a mirror plane perpen- 
dicular to St, will relate: 



and 6«±3" 1 . 

For (uvw)lSt therefore q = q e if G d =2 and y 
if Cfi = m or mm 2. 

In summary, therefore, 

I- N < i 

J -¥ (i - q( '~ l 



(5) 



where g P ^ g. 

For noncubic G* and \S t \ > 2. 

If S ^ || ( uvw ) _ 

If S t 1 (uvw) and S, || 2 
If St J- ( uvw ) and St not 1 1 2 
If 5, A <atw) # wtt/2 



<?e = 2 
Qe = q 
q e = q. 



The value of </ P can be expressed in words as the 
following count taken from the q nonself-inverse 
operator pairs: A pair is counted if neither of its 
elements multiplied by any element of the defect 



symmetry equals either its own inverse or an element 
of a pair previously counted. 

4. Jump Rates for Special Defect Directions 
in Noncubic Point Groups of the Trap Site 

For special defect directions (uvw) in noncubic trap 
point groups G t an exhaustive list of the total number 
J, of distinct paths, is given in tables 2 and 3 relating 
to centrosymmetric and other G t , respectively. The 
maximum and minimum number of jump rates needed 
for relaxation is also given in these tables. If the con- 
cepts involved are not adequately clarified in the intro- 
duction of this paper, the reader is referred to our 
previous publication [1]. It is noteworthy that for 
special positions in noncentrosymmetric point groups, 
«/min>2 occurs only for G ( = 6m2 and Gr/ = m[ ooi]. 

The last two columns of table 2 show additional 
figures in parentheses corresponding to relaxation 
jump rates needed after the population has been upset 
by an influence that itself was centrosymmetric such 
as homogeneous mechanical strain. For general posi- 
tions these numbers are always one less than the corre- 
sponding one for noncentrosymmetric influences as 
previously discussed [1], except for the minimum 
number of paths in point groups having a 3 axis, for 
which the numbers in parentheses did not differ from 



Table 2. Paths for centrosymmetric noncubic trap symmetries; all special positions of point defects 
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'Fable 3. Paths for noncentrosymmetric noncubic trap symmetries; 
all special positions of point defects 
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that corresponding to noncentrosymmetric influences. 
For special positions, however, there are additional 
causes for equalities of numbers in and out of paren- 
thesis while in point groups containing a 3 axis a 
difference in minimum numbers of paths is obtained 
when the special direction is parallel to that 3-axis. 
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